Area evolution, bulk viscosity and entropy principles for dynamical horizons 
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We derive from Einstein equation an evolution law for the area of a trapping or dynamical horizon. 
The solutions to this differential equation show a causal behavior. Moreover, in a viscous fluid 
analogy, the equation can be interpreted as an energy balance law, yielding to a positive bulk 
viscosity. These two features contrast with the event horizon case, where the non-causal evolution 
of the area and the negative bulk viscosity require teleological boundary conditions. This reflects 
the local character of trapping horizons as opposed to event horizons. Interpreting the area as the 
entropy, we propose to use an area/entropy evolution principle to select a unique dynamical horizon 
and time slicing in the Cauchy evolution of an initial marginally trapped surface. 

PACS numbers: 04.70.Bw, 04.20.-q, 04.70.-s 



Introduction. A new paradigm has recently emerged in 
the theoretical approaches to black holes, following the 
introduction of future outer trapping horizons (FOTH) 
by Hayward and that of dynamical horizons (DH) 

by Ashtekar and Krishnan Q (see Refs. 0,0 for reviews). 
This new approach relies on local characterizations of 
black holes, via trapped surfaces, and contrasts with the 
traditional treatment which is based on the notion of 
event horizon (EH). The latter is a highly non-local con- 
cept and requires the knowledge of the whole spacetime 
to be determined. This feature makes the EH a not very 
practical representation of black holes for studies in nu- 
merical relativity and quantum gravity, and motivated 
the new approach. 

In this note, we focus on the evolution of the area of 
cross sections of FOTHs and DHs. We obtain an area 
law of different nature from the "area increase law" de- 
rived by Ashtekar and Krishnan 0, • Inspired by the 
membrane paradigm developed for EHs 0, > we treat 
the cross sections of FOTHs and DHs as viscous fluid 
bubbles and interpret the area law as an internal energy 
balance equation, completing the Navier-Stokes-like mo- 
mentum law derived previously [jj. It is then apparent 
that the bulk viscosity is positive. Associating area with 
entropy, we also discuss the use of an entropy principle 
to pick up a unique dynamical horizon in an evolution. 
Both a maximum entropy production and a time con- 
vexity requirement (interpretable in terms of Clausius- 
Duhem inequality) are considered. 

Evolution of the horizon area. Both the EH and the 
FOTH/DH approaches to black holes can be formu- 
lated in terms of a hypersurface Tt embedded in a four- 
dimensional spacetime (M.,g), and foliated by a family 
of closed (topologically S 2 ) spacelike 2-surfaces (S t )tm- 
If it represents an EH, TL is a null hypersurface and any 
foliation (Stjtes. is admissible. In the case of a DH, TL is a 



spacelike hypersurface and the foliation (iS t ) tg R is unique 
[t defined up to a relabeling t \— ► t' — f(t)], each St being 
a marginally trapped surface 0. 

Given the foliation (St)teM of W, there is a unique 
"time evolution" vector h that is tangent to TL, orthogo- 
nal to iS t and obeys Ch t = 1, where Ch denotes the Lie 
derivative along h: Cut — h^d^t. The latter property 
implies that the 2-surfaces St are Lie dragged by h. Let 
C be half the scalar square of h with respect to the met- 
ric g: h ■ h = 2C. It is easy to see that the sign of C gives 
the signature of the hypersurface TL: C is positive, zero 
and negative for respectively spacelike, null and timclike 
hypersurfaces. There exists a unique pair (£., k) of null 
vectors normal to St and a unique vector m normal to 
TL such that 



h = £-Ck. 



£ + Ck and Ik 



-1. (1) 



For any vector field v normal to St, such as h, m, t or fc, 
we define the expansion 9^ and the shear tensor o-W 
of the surface St when Lie-dragged along v by 



C v q 



) q + 2c^( ,, ) and tr<r ( ^=0, 



(2) 
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where q is the induced metric on St (q is positive defi- 
nite since St is assumed to be spacelike), C v q is its Lie 
derivative resulting from the dragging of the surface St 
along the normal vector v, and tr cr( v > is the trace of cr^ v ' 
with respect to the metric q. A consequence of (01 is that 
Q( v ) governs the variation of the surface element 2-form 
e s of S t according to C v e s = 8^e s . 

Let us recall that St is called a trapped surface if 9^ < 
and 6^' < 0, and marginally trapped surface (MTS) if 
6>( fc ) < and 6>W = 0. The hypersurface H = UteR 5 * 
is said to be a future outer trapping horizon (FOTH) 
if (i) each St is a MTS and (ii) the (local) outermost 
condition CkO^ < is satisfied lj. H is said to be a 
dynamical horizon (DH) if (i) each St is a MTS and (ii) 
TL is a spacelike hypersurface 0. Note that in generic 
dynamical situations, the notions of FOTH and DH are 
equivalent [j| . In stationary situations, a FOTH becomes 
a null hypersurface, whereas a DH (which by definition 
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is spacelike) cannot exist; it should be replaced by the 
notion of isolated horizon (IE) 0, IE El El- 

If 7i is an EH, the 2-surfaces St are not MTS, except in 
stationary configurations (Kerr black hole). On the con- 
trary they are expanding, by the famous Hawking area 
increase law 0: 0^ > 0. 

Let us denote by k the component along £ of the "ac- 
celeration" of h in the decomposition [l3| 

V h h = n£+{Cn- C h C)k-VC, (3) 



where V is the spacetime connection and X> the con- 
nection associated with the metric q in St. Besides let 
fl^ be the normal fundamental form of the surface St 
(also called connection on the normal bundle) defined by 
fiW := — k ■ V^£, where q denotes the orthogonal pro- 
jector on the surface St (see e.g. Sec. III.D of Ref. y|)- 

From the Einstein equation, one can derive the follow- 
ing evolution law for any foliated hypersurface Tt (details 
are provided in Ref. 



C h 9 {m) = K (h) - -9^9 {m) -a {h) : <r {m) -8irT(m,h) +0 {k) C h C + T> ■ (2Cn (£) - VCj , (4) 

where T is the energy-momentum tensor of matter (if any), an upper arrow indicates index raising with the metric q 
and the notation ":" stands for the double contraction, i.e. cr( h ' : <j l<rn ^ := a^ m ^ ab . If we specialize Eq. (0} to the 
cases of (i) an EH and (ii) a FOTH or a DH, we obtain respectively 

Ci9& + (0 W ) 2 -KB® = \{9^) 2 - „W : ff W - 8irT(£,£) (5) 
C h 9 {h) + {9 {h) ) 2 + k9 w = ^{9 {h) f + (T {h) : cr (m) + 8nT(m,h) -9 {k) C h C + T>- (:DC* - 2Cft w ) . (6) 

I 



For the EH, we have used the null character of H, which 
implies C = and h = m = £, yielding Eq. JSJ. This is 
nothing but the null Raychaudhuri equation for a surface- 
orthogonal congruence In this case Eq. J3J reduces 
to = k£, i.e. k is the non-affinity coefficient of 

£ and coincides with the so-called surface gravity for a 
stationary horizon, provided that £ is normalized in terms 
of the stationarity Killing vector. For the FOTH/DH 
case [Eq. ©], we have used the property 0w = to write 
0(m) = -ff(h) ^ since E q . gives 6K m ) = -9^ + 20^. 

The area of the 2-surface St is A(t) = J s e s . Since 
the surfaces St are Lie dragged by the vector h associ- 
ated with t, the first and second derivative of the area 
with respect to t are dA/dt = f s 9^e s and d 2 A/dt 2 = 

J St [C h {K) + (0 (h) ) 2 ] e s . Then, assuming C > in 
Eq. 10 (i.e. considering a DH only), we may introduce 
k' := -C~H ■ V h h = k- C h \nC [cf. Eq. ©], and 
integrate Eqs. JSJ and |JBJ over St, noticing that the in- 
tegral of the divergence term in Eq. © vanishes, to get 
respectively 



d 2 A 



dA 



dt 2 K dt 



s, 



(7) 



d 2 A 
dt 2 



,dA 
~dt 



8irT{m,h) +a- (h -> -.a- 



+ ^ Y ^ + (K>- K ')9W]e s , (8) 



where n and k' denote the mean value over St of k and 
k ': R — R(t) := A^ 1 J s ke s (idem for k'). Assume for a 
moment that R and R' are constant and positive, as for 
the EH of a Kerr black hole with foliations compatible 
with the stationarity Killing vector. Let us consider first 
the evolution of an EH, i.e. the differential equation JJJ. 
The general solution of the corresponding homogeneous 
equation is A(t) = consti • exp(St) + const2. Thus, if one 
were solving Eq. as a Cauchy problem, one would ob- 
tain exponentially diverging solutions (since R > 0). It is 
well known that the correct treatment must be teleolog- 
ical (|E 13 or Sec. VI. C. 6 of Ref. 0), i.e. one imposes 
the boundary condition dA/dt = at t — +oo, to get 



dA 
~dt 



D{u)e R{t ~ u) du, 



(9) 



where D = D(t) stands for the integral in the right-hand 
side of Eq. 0. On the contrary, the differential equa- 
tion JSJ for the DH area is such that the solutions of the 
homogeneous equation are decaying exponentially. Ac- 
cordingly, the treatment as a standard Cauchy problem 
from the initial condition dA/dt(0) = Aq leads to the 
non-diverging solution 



dA ■ 
-dt =A ° 



D'(u)e R '<- u - t '> du, 



(10) 



where D' = D'(t) stands for the integral in the right- 
hand side of Eq. (jHJ. The striking difference between 
Eqs. © and (jTU)l is that Eq. ljTU|) is causal (the solution 
at a given instant t depends only on the behavior of the 
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source D' at instants u < t), whereas Eq. © is not. This 
reflects the non-local character of EHs mentioned in the 
Introduction. Of course, in general R and R' are not 
constant (except for small perturbations of a stationary 
black hole [l5| ) , but the behavior described above should 
remain the same. 

Comparing with previous works, one can show that at 
the limit of small departure from an IH, the area law © 
reduces to that established for slowly evolving horizons 
by Booth & Fairhurst 0]. I n the full dynamical regime, 
the evolution law JHJ is different from the area law ob- 
tained by Ashtekar & Krishnan 0,0 (see also ^(|). In- 
deed the latter is derived from a different component of 
Einstein equation: the T(m, I) one instead of T(m, h) 
for Eq. (JHJ). More precisely, the Ashtekar-Krishnan law 
is written in terms of the areal radius R := y A/4ir and 
expresses the variation of R between two surfaces <S tl and 
St 2 (cf. Eq. (3.25) of Ref. 0). Using our notations (cf. 
Table II of Rcf. |S| for the correspondence) its differential 
version (i.e. writing t% = t\ + dt) is 



1 dR 
2~dT 



8n 



s t 



c 



dR 
dt 



1 / „s d R 



(11) 



where n$ := fiW - T>\nC. Notice that this is a first 
order equation in A(t) — 47ri?(i) 2 , whereas Eq. (JSJ) is of 
second order. Moreover, it contains dR/dt on both sides, 
so that one can divide by dR/dt to get an equation which 
does not contain R(t). In this respect, Eq. (|llfl does not 
appear as an evolution equation for A(t). Actually, as 
shown by Hayward , Eq. can be obtained by inte- 
grating over S t the relation expressing that Ch 0^' = 
on TL. This relation involves T(m,£), in the same man- 
ner as Eq. © relates C h 8 {m) and T(m,h). Explicitly 
we have 



,q. (|4|) 



c 



+8vrT(m,^) = 0, 



(12) 

where 1Z denotes the Ricci scalar associated with the met- 
ric q in St- Dividing Eq. (|12fl by 8ttC, integrating over 
St, invoking the Gauss-Bonnet theorem to set the inte- 
gral of 1Z to 8tt and multiplying by dR/dt leads to the 
differential form ljll|l of Ashtekar-Krishnan law. 
Energy dissipation and bulk viscosity. In the membrane 
paradigm approach to black holes, Price and Thorne |7J 
defined the surface energy density of an EH as £ := 
—0^/8it and interpreted Eq. as an energy balance 
law, with heat production resulting from viscous stresses. 
By analogy, let us define the surface energy density of 
a FOTH/DH as e := -6^ /8tt, where the role of the 
normal to TL is now taken by m instead of £. Since 
0(m) = _ (h) for a FOTH/DH, we have e = 0^/8tt 
and we may rewrite Eq. (j^J) as 



C h e+9We 



-— 0( fc > 

8tt 



rim) 



87T 



16tt 



-V-Q + T(m,h)-—C h C, (13) 

57T 



= —^vb, where ~uj 



with Q := 3L COW -1/2'DC 

is the anholonomicity 1-form (or iwist 1-form) of the 
2-surface <S t (see also Sec. IV. A of Ref. ||) and 
denotes its vector dual. It is worth to write aside Eq. 1|13|) 
the generalized Damour-Navier-Stokes equation 



-T> ( — 

\8tt 



"8^T 



-T(m,9) + — ©C, 

07T 



T>Q( h ) 
16tt 

(14) 



where 7r := — Jl £ /87r. It is striking that Eqs. I|13|l - 
(|14|l are fully analogous to the equations that govern a 
two-dimensional non-relativistic fluid of internal energy 
density e, momentum density 7r, pressure k/8tt, shear 
stress tensor cr^ m ^ /8n, bulk viscosity £ = 1/167T, shear 
strain tensor <rW, expansion 6*^-', subject to the exter- 
nal force density —T(m,q) + 8^ /8ir'DC, external en- 
ergy production rate T(mJi) — 0( k > /8w C and heat 
flux Q (see e.g. Ref. 17]). Let us notice that the 
shear viscosity fi does not appear in Eqs. H13|) - (|14fl . be- 
cause the standard Newtonian-fluid relation between the 
shear stress tensor er( m ) /8n and the shear strain ten- 
sor cr( h \ namely <r^ m ' > /8tt — 2ficr^ h \ does not hold. 
Here we have ct^/Stt = [cr^ + 2Cct^]/8tt, so that 
the Newtonian-fluid assumption is fulfilled only if C = 
(IH limit). On the contrary, it appears from Eqs. 113|) - 
(|14fl that the trace part of the viscous stress tensor S v i sc 
obeys the Newtonian-fluid law, being proportional to the 
trace part of the strain tensor (i.e. the expansion 
trS visc = 3C# (h) . 

Let us point out two differences with the EH case |a,l3. 
First the heat flux Q is not vanishing for a FOTH/DH, 
whereas it was zero for an EH. Notice that Q is a vector 
tangent to St so that the integration of Eq. H13|) over the 
closed surface St to get a global internal energy balance 
law would not contain any net heat flux. The second ma- 
jor difference is that the bulk viscosity £ is positive, being 
equal to 1/167T, whereas it was found to be negative, be- 
ing equal to — 1/16tt, for EHs This negative value, 
which would yield to a dilation or contraction instability 
in an ordinary fluid, is in agreement with the tendency of 
a null hypersurface to continually contract or expand, the 
EH being stabilized by the teleological condition impos- 
ing its expansion to vanish in the far future. The positive 
value of the bulk viscosity found here shows that FOTHs 
and DHs behave as "ordinary" physical objects and is in 
perfect agreement with their local nature. 
Entropy evolution and choice of TL. Hitherto we have 
discussed the horizon TL as a given hypersurface in M.. 
We adopt now a 3+1 approach and consider, in the line of 
Ref. [18| , the world-tube evolution of a MTS So contained 
in an initial Cauchy slice So of a 3+1 foliation (Et)teR- 
Denoting by s the unit-normal vector to Sq contained in 
Sq and by n the unit-normal vector to Sq, we decompose 
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H's "time evolution" vector as h = Nn + bs, where N is 
the lapse function associated with the time parameter t of 
the 3+1 foliation. Introducing the null vectors £ := n + s 
and k:=(n- s)/2, we get h = (b + N)/2 1+ (b - N)k 
and C = (b 2 — N 2 )/2. Combined uniqueness results from 
Refs. show that: (i) to each 3+1 foliation (E t ) tg R 

there corresponds a unique DH containing So and sliced 
by MTSs S t C and (ii) different 3+1 slicings lead 
generically to different DH. In other words, the evolution 
of Sq into MTSs is an ill-defined concept, since different 
DHs pass through Sq. A natural question consists in 
introducing a notion of preferred DH or, equivalently, of 
a preferred 3+1 slicing. 

With this ultimate aim, let us push forward the viscous 
fluid interpretation of FOTH/DHs. As an ever growing 
quantity on TL, it is natural to interpret the area A as 
proportional to the entropy of the horizon, the equilib- 
rium case providing the relation S = ^-A 0,0, 151. In 
this sense, Eq. (JSJ) for an EH has been interpreted in the 
membrane paradigm 0] as a viscous fluid entropy dis- 
sipation equation for the differential element of entropy, 
with the specificity of involving a second time derivative 
term in the entropy production. Instead of the analogous 
discussion for FOTH/DHs based on Eq. ©, we take the 
condition C^G"' = 0, i.e. Eq. (|12|) . as a starting point. 
When expressed in terms of 3+1 quantities, it becomes 
an elliptic equation on (6 — TV) with a source proportional 
to (b + N) [lillijj- In this 3+1 approach, a second rela- 
tion between (b — N) and (b + N) is necessary in order 
to determine 7i. 

In a first attempt to get this second relation, we con- 
sider a maximum entropy production criterion, which can 
be motivated in the context of non-equilibrium thermo- 
dynamics . Maximizing dS/dt leads to [l3l ] 

b- N= -const- 6^. (15) 
In loose terms, and independently of thermodynamic con- 



siderations, Eq. (|15fl singles out the DH that approaches 
"the fastest" to the EH. However it leads to a function 
A{t) which is only C° in the matching with an initial IH. 

Alternatively, we can control the response of 7i to the 
arrival of energy/matter, by prescribing the convexity of 
the entropy in time. We propose a phenomenological 
choice for the second derivative of the area element a := 
C h 6^ + (9<-V) 2 : 

0(h) 0(h) 

a = F[sources] + a—^-C h C + f3—^-C h C h C, (16) 

such that F [sources] > whenever matter or gravita- 
tional radiation crosses the horizon [e.g. T (£,£.) + 1 
or er^ =/= 0], and F[sources] < otherwise. Inserting 
Eq. (|16|) into Eq. JHJl leads to a second-order evolution 
equation for C = (b 2 — N 2 )/2, whose resolution provides 
the additional relation between (6 — TV) and (b + N). Ini- 
tial condition C = guarantees the C 1 matching with 
an IH, whereas the choice of initial Ch C amounts to the 
choice of DH (consistency requires an appropriate choice 
of parameters a, (3). Alternatively, the equation for C 
can be seen as a (now first-order) balance equation for the 
element of entropy. Interpreting it as a non-equilibrium 
thermodynamics Clausius-Duhem-like inequality pol ] (by 
enforcing the positivity of the entropy production), pro- 
vides a guideline for fixing the phenomenological term 
F [sources]. Other phenomenological options to Eq. ljTrj|) 
can be proposed, and we aim at exploring them numer- 
ically [I3 A tempting possibility is to base the (un- 
avoidable) choice of DH, upon an entropy principle de- 
rived solely from the structure of the hyperbolic system 
defined by (part of ) the Einstein equations on TL. 
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